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Abstract: The Chapman-Richards Function and its two exception cases in applications were discussed and compared with 
the Schnute model in stand growth studies. Compared from ail perspective, it was found that the Schnute model commonly 
used in forestry was identical to the Chapman-Richards function. If some parameter in the Chapman-Richards Function was 
unconstraint, the function could also be very versatile to tit some exceptional growth curves, the fitted function should be 
identical to that the Schnute model. 
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Introduction 

As a growth' model for an individual tree or forest stand, 
the Chapman-Richards Function , an extension to the 
Von Bertalanffy’s growth theory( Bertalanffy 1957) 
embodies all other established growth function. This 
function has been widely used in forestry because of its 
flexibility, accuracy, and meaningful analytical proper¬ 
ties (Cooper 1996, Pienaar and Turnbull 1973, Osumi 
and Ishikawa 1983, Ito and Osumi 1984) 

Based on a simple biological theory of accelerating 
growth, Schnute(1981) developed a new and compre¬ 
hensive growth model with statistically stable parame¬ 
ter, that was said to cover all other growth functions 
existing to date. The Schnute model seems to have been 
accepted quickly because of its versatility and stable 
parameter estimation(Bredenkamp and Gregoire 1988, 
Yang and Feng 1989, Zeide 1993). Bredenkamp and 
Gregoire) 1988) and Yang and Feng(1989) studied the 
performance of the Chapman-Richards Function along 
with the Schnute model. They concluded that the com¬ 
petition mortality, the stand resumed growth beyond 
what appeared to be an asymptotic level where the 
Chapman-Richards Function did not fit the data satis¬ 
factorily, Meanwhile, the Schnute model can tract such 
a trend in growth successfully. 

The above conclusion was drawn probably because 
the model was only used for data fitting, biological and 
local asymptotic properties of parameters might have 
been neglected. In this paper we demonstrate that the 


Chapman-Richards function and the Schnute model 
have the same capability, they can fit the same data set 
equally well and with the same accuracy. The objective 
of this study is to compare the two models to show the 
relation between them. 

The Chapman-Richards Function 

In studying plant growth, Richards) 1959) proposed to 
extend the range of parameter m in Von Bertalanffy’s 
growth model to m > 0 rather than m >1 , i.e., 

^ = rp/"-K y (]) 

dr 

Where y is a measure of plant volume, t is plant age, q, 
k and m are parameters. It is also assumed that q and k 
<0 if m >1, q and k > 0 if 0 < m <1. By integrating 
equation (1) with the initial condition of y=y 0 for t=0, 
the Chapman-Richards function is obtained as 

I 

y = A(\~Be- a y-, 0<m< 1 

, ( 2 ) 

y = A(] +Be m> 1 

with 

B i », 

n 1 “ ^ 

/f = (-)' B = ± - , K = (]-m)K 

k q 

K 

where A -is tjpe asymptote, m determines the shape and 
inflection point of the growth function, and K is related 
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to growth rate, equation (2) has the following proper¬ 
ties: 

• It has asymptote A. y->A as t-»3c. 

• y is a monatomic increasing function of/ 

• The inflection point is at: 

' B 1 

/. - ln( -). v, = A m 1 when ()</;/< I (3) 

' K I - m ’ 


Equation (7). a species case of the Chapman-Richards 
function, can describe the growth pattern of a stand 
resuming growth beyond what had appeared to be an 
asymptotic level after marked competition mortal- 
ityfBredenkam and Gregoire 1988). 


/, = - ln( v, = A m 1 when m> 1 

' K I - m ' ' 

• The curve has a cross-over point at: 

/,, = ln(/i) (4) 

A 

such that t —0 at /=/„ 

In addition to m > 0. the following two exceptions 
could occur in actual applications. 

Case 1 

for <0 and q. k >0. equation (2) was fitted to some 
data of fast-growing species with no inflection 
point, fhe trees continue growing from an early age as 
shown in figure 1 and the corresponding growth func¬ 
tion should be: 

I 

i .1(1 Be kl )' with A, B. A > 0.///<() (5) 

Equation (5) does not have an inflection point because 
m<0(see y, in equation (3)). 

Case 2: 

for in <0. rpO. and k< 0. let k'=-k. than (1) become 


The Schnute Growth Model 


fhe Schnute model used in forestry is developed as 
follow: Let )’(/) be the size of an organism with growth 

rate d)7d/ and relative growth rate /. = ■■■ ^ . The 

)' Ji 

relative growth rate of Z. . is assumed to be a 

Z di 

linear function of Z as: 

^ - -(a t b'/.) (10) 

Z d / 

where a * 0 and h * 0 are parameters. Since a *0 and h 
* 0. the growth could, be accelerating or both. (10) has 
another form: 


(-.- - )d/ • a d/ (II) 

a s hZ Z 


Integrating (II) with the initial condition of / = l\ and 
Z(T | )=Z| is an unknown constant gives: 


, a + hZ Z. 

In (-■-- / ) = a(l-T,) 

Z a + hZ, 


from (12) Z(/) can be written as: 


( 12 ) 


: = rp-"' + k‘ 1 v. q.K ’> 0. /// < 0 (6) 

d/ 

fhe solution to (6) with initial condition is: 

I 

v A’( Be K ' - I) 1 A'. B', C> 0 ,m <0 (7) 

w here 

A 1 - ( n )' B'= 1 +('--)' A' 1 = k’(1 -/;?) 

k’ A' 

Equation does not have an asymptote. After the curve 
passes the point (/,. v, ). dt/d/ begins to increase as t 
increase. Such a growth curve is analogous to the one 
shown in figure 2. The inflection and cross-over points 
for (7) are: 

1 I fll ' 

i, ----( )ln( m ). v, =A'( 1 -m)' (8) 

A' B' ’ 


d ln( Y) 
di 


/ ./(/ 
C//.C’ 

= /(/) = - 

a + hZ . (1 -c 


i a 
no 


(13) 


I d 
b d / 


ln(u +/?Z,( I - e 


The solution to (13) with t-r, at t=7', is: 


a+hZA \-c "" 11 

) (o v,(- )'' 


(14) 


a 

where integration contant Z, can be solved using the 
terminal condition of / = 7\ and )'(t)= is as: 


= 


ci(yS — .v,'') 
h\\ h ( I - e ' ! ') 


(15) 


Substituting Z, into (14) gives the explicit solutions to 
(10) or the Schnute model as 
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1 _ ? olt * 1 
Y(n = (y:+(y 2 h -y:) r ~-^ r -,~V 
I — e 


(16) 


where T ] and T 2 are initial and terminal ages, y, and y 2 
are corresponding size at 7) and 74. 

By definition. /„ is the value oft such that )'(/)=0. Let 
(16) equal to 0. is derived as: 


- 7; + 7t 


I c‘" : v/ ~e ,,r '\\ 

-ln ( - -T—- -1- 

u w - y. 


(16) has the upper asymptote as: 

a/, i' ,ii. 1' 1 

e — e y, h 

3 ' ~ I ~ ..37: I 

c' - e 

The inflection point is at: 


(17) 


(18) 


Z = c,e ’"(a + hZ) (24) 

and 


y = c ; (l - c\be ‘") h (25) 

where c, and c 2 are integration constants. From equa¬ 
tion (24) we can get c x e Z / (a + bZ ), substituting 
this into (25) generates: 


h b U 
= C\ - 

a + hZ 

Since Z=( l/y)(dy/d/), (26) can be written as: 

dr a hit, u 
— = — c\ v — v 

d / h h ' 


(26) 


(27) 


* / / til. b al\ h x 

... T I . Me -y. -e y r ) 

7 ‘ 7. hit .— )• 


v. 


6)(e__v ; - • e' v ) > 


Let m = \- h. K~a/h,x\=(alb)c 2 h then (27) becomes (20), 
which is the differential equation for the Chapman- 
(19) Richards function. The above exercise proves that the 
two model have same differential equation with uncon¬ 
strained parameters. 


Relationships between the Two Models 

By comparing assumptions, solutions and relations 
among the parameters for the two models, one can 
show that they are interchangeable. 

Identical differential equations 

If q. k. and m in equation (1) are unconstrained or: 

- i\v ni -Kv. q. k. m are anv constant (20) 
d 1 

multiplying both sides of (20) by y' 1 gives: 

I c l'’ 1 

Z(/)= -'-=iy -k (21) 

y di 

Differentiating (21) by d l gives: 

^ = q(m - Dv" : = (in - Dili’"' '7 (22) 

d / d / 

liquation (21) indicates that r\v"''' = Z + k. substitut¬ 
ing this into (22) gives: 

= -(k(1-;h) + (1-/;i)Z) (23) 

Z d 1 

Let a=K7l-/»; and 6=1 - 111 . then (23) becomes (10) or 
the differential equation for the Schnute model. The 
Schnute model can be derived from the assumption for 
the Chapman-Richards function and vice-versa. From 
(132) and (B5) (see Appendix B). there exist: 


Identical solution set for the general solutions 

It is shown in Appendices A and B that the general 
solutions to the two differential equations have the 
same formulation, their parameters are interchangeable. 
For example, the parameters in the general solution for 
the Chapman-Richards function (A5) cambe written as: 

/n=l-/),K=-,q= (—)c,, c, = ac t c , 
b h 

where a. b ,and c 2 are parameters in the solution for 
the Schnute model (B5). Similarly, a. b, c x and c 2 can 
be expressed as: 

I 

a =k(\ - 111 ), h = 1 - m, c\ =(—)'"'”, c. = -- 

i\ q(l - 111 ) 

Identical explicit solution 

The general solution to (20) is: 

I 

v = (- + ce K . ’)'"' (28) 

K 

where c is the integration constant. With the initial con¬ 
dition of 1 = 7’,and Y(t) =y,. c becomes: 

c = (y, 1 " - ’-K" W]L (29) 

K 

Putting c back into (28) to have: 


139 





Comparison of the Chapman-Richards Function with the Schnute Model in Stand growth 


7(f) = (- 1 ! + (y '- m - H)e(30) 

K ' K 

t)/k can be solved with the terminal condition of t = 7V 
and 7(f) = y 2 . Substitute t)/K into (30) and use b= (1-m), 
a= k(1 - m), then the solution to (20) is: 


decided by m and Kora and b. The asymptote exist if 
0</w<l (0<A<1) and k>0 ( a >0) or m <0( b> 1) and k 
>0(a >0). If m > I. the asymptote exists only if 

-/>ln( V -) 

*>- 

t 2 -t, 


W = (y,*+(y2*-y,*) l 1 _*-^Z ) i (3D 

which is the Schnute growth model. 

The above analysis reveals that if the parameters in 
Von Bertalanffy’s model are not constrained and the 
two stages are given for a growth process, then the 
Schnute model can be deduced from (20) by using the 
initial condition of (T h 7,) to get a specific solution first, 
then the terminal condition of (T 2 , y 2 ) is used to replace 
t|/k (in equation (30)) which is related to the asymptotic 
parameter A{x\Ik=A u ") in Chapman-Richards function. 
In this case, the Schnute model is an explicit solution to 
(20) with an implicit asymptotic parameter. 


The relation between the cross-over and inflection 
points 

The cross-over point of the Schnute model can be de¬ 
rived from that of the Chapman-Richards function: 


1 


K a e - v\ -e 'y x 

t <//, b a/ h 

a y 2 — y i 


= T. +T, 


In the same fashion, (19) can be transformed from (3), 
implying points)/,, v,) and (/*, 7*) are identical. The 
relation between asymptotes A and y, is obvious. 


Examples 

Two examples are used to verify the above results. 

Example 1 

Yang and Feng (1989) used both model to study the 
volume growth of permanent sample plots of Cryp- 
tomeria spp with 5 initial densities. Table 1 consists of 
the parameters by Yang and Feng (1989). With the 
given relations between the parameter, parameters A, B, 
k, and m in the Chapman-Richards function could also 
be calculated from the value of y u y 2 , a, and b in the 
Schnute model. The calculated results are the same as 
the fitted ones. Since 0</n<l(or)</h<l), the inflection 
point, asymptote, and /„ exist for all 5 initial densities. 


Density 

A 

B 

K 

m 

The Chapman-Richards function 

RSS DF 

MSE t„ 

/, 

V, 

3.1 v 3.1 

567.5702 

0.8291 

0.1510 

0 8907 

918.643 

23 


39.9410 -1.2413 

13.4197 

197 

2.3 x 2.3 

644 397 

1 4379 

0.1465 

0.7362 

8077.233 

23 


351 1841 2 4791 

11.5750 

202 

1.9 x 1.9 

804.9043 

1.5185 

0.1146 

0.6328 

1079.113 

23 


46 9180 3 6301 

12.3726 

231 

1 7 x 1.7 

796 7438 

1 5093 

0 0865 

04820 

11341.77 

23 


493.1205 4.7608 

10.8283 

175 

1.5 x 15 

1285 9725 

1.2163 

0.0395 

0.0827 

10650.108 

23 


463.0482 49590 

7.1457 

85 

The Schnute model 

Density 

V, 

V, 

VS 

a 

h RSS 


DF 

MSE 1,. 

t* 

• V* 

3.1 x 3.1 

567.5702 

6.1977 

498.57 

0.1510 

0.1094 918.643 


23 

39 9410 -12413 

13 4197 

197 

2.3 x 2.3 

644 369 

7.4926 

579.65 

0.1465 

0.2383 8077.233 


23 

351.1841 2.4791 

11.5750 

202 

19 x 1.9 

804.9043 

4.2045 

622.26 

0.1146 

0.3672 1079 113 


23 

46 9180 3 6301 

12.3762 

231 

1.7 x 1.7 

796.7438 

11152 

609 99 

0.0865 

9 5918 11341.77 


23 

493 1205 4.7608 

10 8283 

175 

1.5 x 1.5 

1285.9725 

1.1668 

711.66 

0.0395 

0.9173 10650.11 


23 

463.0482 4.9590 

7.1457 

85 


Note (1) T1 =7 and T2= 27. (2) t is stand age(vear) and y is stand volume (in') 


The relation between two sets of parameters 

Parameters a and b in the Schnute model can be repre¬ 
sented by those in the Chapman-Richards function as: 

a = K - (1 - in )k 
b = 1 - w(or m = I - b) 

Similarly, there are: 


■4 = (±(.7 + /\-iT„, » 


+7;> 


B = ±- 


1 - <?" 

, h h , 

(7: -Y, ) 


oi. b 

e--y 2 


■ e“'' \\ 


such that A >0 and B >0 

Whether the two model have the asymptote or not is 

Table 1. The parameters estimated for equation (2) and (16) 
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Example 2. 

Bredenkamp and Gregoire(1988) fitted both models to 
the mean diameter growth data of Eucalyptus grandis 
plantations with 12 initial densities. Table 2 contains 
the estimated parameters and mean squared residuals 


(MSE) from both models (Bredenkamp and Gregoire 
1988) with some re-arrangement. Missing values in the 
table indicates the case where the Chapman-Richards 
function failed. 


Table 2 Estimated coefficients from both-parameter Chapman-Richards and Schnute models to describe the dbh development 
of Eucalyptus grandis stand( Missing values indicate stand densities where the Chapman-Richards model could not be adequately 


fitted (Bredenkamp and Grgoire 1989). 

Density 

A 

B 

Chapman-Richards 

K m 

MSI: 

Vi 

V 2 

Schnute 

a 

b 

MSE 

25 

950.56 

1.0632 

0.0938 

0.1948 

2354.61 

38.90 

893.18 

0.0938 

0.8052 

2354.61 

62 

838.15 

1.0989 

0.0883 

0.0204 

208.78 

29.29 

786.38 

0.0883 

0.9796 

208.78 

124 

712.09 

1.1168 

0.1030 

-0.0074 

189.90 

31 43 

685,26 

0.1030 

1,0074 

189.90 

247 

559.43 

1.1365 

0.0973 

-0.3207 

133.09 

26 51 

539,57 

0.0973 

1.3206 

133.09 

371 

489.55 

1.1302 

0.0877 

-0.5913 

218.17 

25.66 

469.80 

0.0877 

1.5916 

218.17 

494 

420.03 

1,1439 

0.0955 

-0 7161 

90 39 

26.53 

407.72 

0.0955 

1.7161 

90,39 

741 

370.84 

1.1543 

0.1014 

-0.8047 

106 10 

26.52 

362.26 

0.1014 

1.8046 

106.10 

988 

338 65 

1.1129 

0.0741 

-1 1164 

168.72 

25 41 

322.57 

0.0741 

2.1159 

168.72 

1482 

326.95 

1.0556 

0.0368 

-1 7056 

96.05 

27.25 

284.30 

0.0368 

2.7056 

96.05 

2965- 






32.96 

272.73 

-0.0247 

3.4050 

99.64 

4304 






30 49 

258.63 

-0.0590 

3 6780 

71.75 

6726 






29 90 

284.24 

-0.1076 

4 1654 

71,76 


Note: Bredenkamp and Gregoire (1989) used c„. c, c 2 . and c, as the parameters for the Chapman-Richards function, in Table 2. A=c„, 


K=c | m=l-l/c 2 and R=c y The DBH was measured in mm 


The estimated parameters for the two model match 
and the MSE's are the same for the fist 9 densities. 
There are inflection point, asymptote, and t 0 for the first 
two densities since 0<m<l(0</t<l). As for densities 3 
to 9 (124-1482 stems /hrrf), only asymptote and t 0 exist 
because m < 0 (/><1), which is the first exception for the 
Chapman-Richards function or (5) (see Fig. 1). 


fitted(Bredenkamp and Gregoire 1988). This situation 
is graphically portrayed by Fig. 2. For the Schnute 
model, there are b>\ and A< 0. It seems the asymptote 
and inflection point do not exist, only tO could be found. 
As point out by Schnute(1981), his model may not have 
inflection point for b> 1, however, the inflection point 
does exist if: 



Fig. I. The growth curve predicted by equation (5). 

The parameters used are those for the density of 124 stems/hnr in table 2 


The Chapman-Richards function could not fit the data 
for the last three densities (2965-6726 stems/ hnf ) 
because after competition mortality, the stand resumed 
growth and the model could no longer be satisfactorily 


-Mn^ 2 ) 

- ^<a< 0 

t 2 -t x 

Parameter a for the last three densities was tested to 
satisfy the above condition. For example, the value of a 
for the density of 6726 stems/hm 2 is 


284 ">4 
-4.165 ln( ' ) 

_ 29,90 

32.83- 1.50 


-0.2994 <a - -0.1076 < 0 


with T ] =1.50 and T 2 =32.83(Bredenkamp and Gregoire 
1988). Therefor, not only the curve / 0 =1.48, it also has 
the inflection point at: 


1 .. b ul\ b 

,*= T t +T, - — ln( ) = 14.74 


a 


y 2 -y i 
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Fig. 2 The growth curve predicted by (7). 

I he parameters used are those for the density of 2%5 stems/hnr 
in fable 3 

It is determined that the following is common for the 
last three densities, i.e.. for the Sclmute model, there are 
a <0 and h> I. for the Chapman-Richards function. m<0 


and k>0. Equation (2) does not describe such data very 
well, the right choice for these stand should be (7). 
Taking the density of 6726 stems/hnr as an example, 
by the relation between the parameters, the Chpman- 
Richards function should be: 

! 

}’(/)= 127.5245 (0.85297c“ ln7 '‘' - l) 1 " JU,M| 

This model has /„ and an inflection point, but no as¬ 
ymptote. From (8) and (9), it is known that / f ,= 1.48 and 
f,=14.74, these results coincide with that from the 
Sehnute model. 

Table 3 contains the values for parameter A , B. k. 
and m for the three densities, these values are calcu¬ 
lated from y,. is. a, and h for the same densities. With 
such estimated parameters. (7). the second exceptional 
case of the Chapman-Richards function, is used to pre¬ 
dict the diameter for the three densities along with (16) 
or the Sehnute model, the predicted values are listed in 
Table 4, showing the two produced identical results. 


Table 3. The value parameter A, B, k, and m calculated from p,, p,, a and b for the last three densities 


* densit) 

.1 

Chapman-Richards 

H . K 

m 

V| 

\\ 

Sehnute 

£/ 

h 

MSI. 

2965 

2604.3V 

04)645 

-0.0247 

-2 4050 

32.96 

272 1a 

-0.0247 

3 4050 

00 64 

4304 

163.X 167 

04)172 

-0.0590 

-2.67X0 

30 4 V 

258 63 

-0.0590 

3.6780 

71 75 

6726 

127 3245 

0.8530 

-0.1076 

-3.1654 

204)0 

2X4 24 

-0.1076 

4 1654 

71.76 


Table 4. the predicted value for the last three densities equa¬ 


tion (7) and (16) 


IX'llsitN 

age(\ear) 

2965 su 

(7) 

rms/hnr 

(16) 

4304 stems/hirv 

(7) (16) 

6726 stems/hnr 

(7) (16) 

2.0 

73.26 

73.30 

64.20 

64.30 

64.32 

64.33 

3 0 

100.18 

100.20 

86.33 

86.33 

84.26 

84 26 

4.0 

116.51 

1 16.52 

99.76 

00.76 

06 40 

96.40 

S 0 

128 03 

128 04 

i lo.io 

1 10.10 

105 88 

105 88 

6.0 

139.23 

130.24 

118 80 

1 18.81 

113 00 

1 14.00 

7.0 

148 16 

148 17 

126 48 

126 48 

121 20 

121.29 

8 0 

156 15 

1 56 16 

133.45 

133.45 

128.05 

128 05 

9.0 

163 42 

163.43 

130 01 

13001 

134 45 

1 34.45 

10.0 

170.15 

170.16 

14 4 00 

145 00 

140 60 

140.60 

1 I o 

1 76.44 

176 45 

151 77 

151.78 

146 57 

146.58 

12.0 

1 82 37 

182.38 

157.33 

157.33 

152 43 

152.43 

13 0 

188 01 

188.01 

162 60 

162 60 

15X21 

158.21 

14.0 

193.39 

103 30 

169 91 

167.01 

163.04 

163 94 

15.0 

198 55 

1 OX 56 

173.01 

173 01 

169.66 

169 66 

16.0 

203.53 

203.54 

17801 

178 01 

175 38 

175.38 

17 0 

208.35 

208.35 

! 82 03 

182.03 

ISI 13 

181 13 

1S.0 

213.02 

213.03 

187.70 

1X7.70 

186.02 

186.92 

19.0 

217 57 

217.58 

102.61 

192.61 

102 76 

192.76 

20.0 

222 01 

222.01 

107.38 

107 38 

198.66 

198.66 


In summary, the Chapman-Richards function will fit 
the last three densities, and the fitted results should be 
identical to those from the Sehnute model. E quation (3) 
and (7). tw'o exceptional cases of Chapman-Richards 
function, are capable of tracking the growth after the 
mortality induced by competition. 

Conclusions 

The Sehnute model with a ^0 and b^O may be regarded 
as an existing model since its assumption, the form of 
solutions, and the cross-over and inflection points with 
which determine the curve shape are identical to those 
for the Chapman-Richards function. In fact, the 
Sehnute model is one of the solutions to (I) if the latter 
is used to model a growth process with two known 
stages. 

Compared with the Sehnute model, the Chapman- 
Richards function is simple and convenient in both 
formulation and its cross-over and inflection point func¬ 
tions. The parameters in Chapman-Richards function 
also have a clear biological interpretation. 

When parameters in (I) are unconstrained, the Chap¬ 
man-Richards function can also be very versatile. Any 
data set which can be fitted bv the Sehnute model 
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should also successfully fitted by the Chapman- 
Richards function. 
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Appendices 


(A) The general solution to the Chapman-Richards 
function 

I he differential equation is 
dr 

; uf -so (AD 

d/ 

Multiply both sides of (A I) by (1 -/;/)t it becomes: 

(I iii)y 11,1 (I - /;/)(r| — kt "') (A2) 

d / 

I et v r 1 '" 1 . then there is 

d.v dr 

(I >n)v 

d/ il / 

Substitute the above (A2) gives: 

d.v d.v 

(1 m){\] K'.v) or (I ;;/)d/ (A3) 

lit I] KV 

I he solution to (A.3) is 

Inti) kAt) ■ k(I m)l + c (A4) 


where c is intergration eonatant 

Placing, v i' 1 '"’ into (A4). the general solution is: 


'1 _ C 

K ki 


.. I I ’’ , K I I "1)1 \ | 

\ ( “ <3 ) . t , = (' 


(A5) 


I et 1 1 1 . [I 


. y k (I - m). o I . can 

k 1 - III 


he simplified as: r - (</ + |tc ' )" 


(B) The Cieneral solution to the Sehnute model 

The differential equation is: 


1 

= -(a + h/) 

/ dt 


<B 1 ) 


The general solution to (HI) is: 


= c 1 . C is integration constant 

a + hZ 


Since / 


1 dr 
r d / 


. (B 2 ) can he written as: 


<B2) 


I dr orp' 

r d/ I hc,L‘ 

(B.3) can be integrated as: 

ln( r) = ' Inf 1 -/vy )tc" 
h 

where c" is integration constant. 
The general solution to (BI) is: 


(B.3) 


(B4) 


r ; c\( 1 -hc t c ) h . r ; = e‘ (BS) 

l et ex =(.■/ . f! = -he, l \ h . y =-a. d = I / b. (BS) can 
also be simplified as r - (a + fir ')" 
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